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. Abstract 

QFT with an external background can be considered as a consistent model only 
if backreaction is relatively small with respect to the background. To find the corre- 
sponding consistency restrictions on an external electric field and its duration in QED 
Q ■ and QCD, we analyze the mean-energy density of quantized fields for an arbitrary con- 

Q_) ' stant electric field E, acting during a large but finite time T. Using the corresponding 

, asymptotics with respect to the dimensionless parameter eET 2 , one can see that the 

leading contributions to the energy are due to the creation of paticles by the electric 
field. Assuming that these contributions are small in comparison with the energy den- 
sity of the electric background, we establish the above-mentioned restrictions, which 
determine, in fact, the time scales from above of depletion of an electric field due to 
, the backreaction. 

■ PACS numbers: 12.20.Ds,11.15.Tk,11.10.Wx 
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I. It is well-known that QFT in an external background provides an efficient model for 
the study of quantum processes in those cases when some part of a quantized field is strong 
enough to be treated as a classical one. For example, QED with an external electromagnetic 
background formally arises from extending the QED Lagrangian by the interaction of the 
matter current with a given external electromagnetic field A^, which is not quantized. 
This is naturally implied as a certain approximation. The study of some problems in QED 
and QCD with superstrong external backgrounds and their applications to astrophysics and 
condensed matter (to graphene physics) has once again raised the question of a consistency of 
such theories. Obviously, the question must be answered, first of all, in the case of a constant 
external field. Calculations that have an immediate relation to the above-mentioned problem 
have first been carried out by Heisenberg and Euler in the case of QED with constant parallel 
electric E and magnetic B fields, see [T]. They computed the change of the vacuum energy 
of spinning particles for arbitrary B and a weak electric field E <C E c = m 2 /e (h = c = 1) 
which is unable to effectively create pairs from vacuum. They interpreted this change as a 
change of the energy of the external field itself and, at the same time, as a change of the 
Maxwell Lagrangian £(°) = (E 2 ~ B 2 ) /8tt by a certain addition C {1 \ The limiting case of 
a strong magnetic field (B ^> m 2 /e, E = 0) yields 



3ir m 

where a — e 2 is the fine-structure constant. This result is in agreement with more advanced 
calculations carried out by Ritus [2J , who arrived at the conclusion that the loop expansion 
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makes sense only for the magnetic fields restricted by the condition B -C F max , 

f max ^exp(^^fO- 
e \ a J e 

Shabad and Usov [3] have recently established a more rigid limitation for the maximal 
admissible strength of magnetic field, B <C 10 28 m 2 /e, having analyzed the structure of QED 
of vacuum in this field, taking into account the interaction of virtual electron-positron pairs. 

The addition has been generalized in a certain way to an arbitrary constant field 
(to arbitrary E) and is now called the Heisenberg-Euler Lagrangian (HEL), for a review, 
see [4]. However, its physical meaning for a strong electric field is not completely clear. The 
problem has been approached from another angle by Schwinger [5], who supposed that the 
variation of an effective Lagrangian of electromagnetic field should arise due to a non-zero 
(in external fields) in — out vacuum current of charged particles, without any restrictions on 
the intensity of electric field; in doing so, he presented its calculation in a constant external 
field, and obtained precisely the HEL. In the general case, the HEL is complex-valued, 
its imaginary part determines the probability of pair-creation, as has been confirmed by 
independent calculations [6l [7]. Considering the in — out vacuum current, Schwinger has 
made it possible to obtain an elegant expression for his effective Lagrangian in terms of 
the causal (Feynman's) Green function. Nevertheless, such an effective Lagrangian is not 
related to the problem of mean values, in particular, it does not reproduce the change of the 
vacuum energy of spinning particles for arbitrary E; the latter problem has to be formulated 
independently as a mean-value problem, and is expressed via noncausal Green's functions; 
see 0. 

In a strong electric field (E m 2 /e, B = 0), the real- valued part of HEL, describing the 
effects of vacuum polarization, is given by the right-hand side of ([T]), where B is replaced 
by E. From this expression one can extract the negative-valued additive correction £w = 
Re to the classical Maxwell density of energy, = D- '. Consequently, in order that 
the total energy density of electric field £ — £^ + £( 1 ' should be zero, the only effects of 
vacuum polarization at the value E ~ F max are sufficient by themselves. On these grounds, 
in [H] it was suggested that F max should be also a limiting value of electric field. It turns out 
that Iml' 1 ' ~ E 2 , which can be interpreted as an evidences that the vacuum instability is 
less than the vacuum polarization. This is not true since the vacuum instability is a nonlocal 
effect, being directly dependent on the electric field duration. 

In our opinion, the most adequate object, whose analysis can answer these questions, 
is the mean value of the energy-momentum tensor of matter, computed with respect to 
various initial states. A detailed calculation of such a mean value in QED in the one-loop 
approximation, taking an exact account of the interaction with the electric background, has 
been given in j9j. Here, we will only use some of these results, in particular, the mean 
energy density of matter for large values of strength and duration of the electric field. We 
examine two cases, when the initial state of the quantized Dirac field is vacuum, and is in 
thermal equilibrium. In addition, we consider the case of charged bosons and QCD with 
an external chromoelectric field. We demonstrate that under these conditions, the effect of 
particle-creation is precisely the main reason for the change of the energy of matter. Making 
a comparison between the change of the energy density of matter and the energy density of 
the external electric field, which is responsible for this change, we obtain restrictions on the 
intensity of the external field and its duration, which we call consistency restrictions. 

II. A constant electric field acting during an infinite time creates an infinite number of 
pairs from vacuum even in a finite volume. This is why we choose the external background 
as a constant electric field acting during a finite period of time; we refer to this field as 
a T-constant field. The finiteness of field duration is a natural regularization in the given 
problem; on the other hand, it is a necessary physical parameter, which subsequently enters 
the consistency restriction on the value of the maximal electric field. The T-constant field 
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turns on at — T/2 = t\ > ti n and turns off at T/2 = ti < t out . We choose the nonzero 
T-constant field potential As(t) as a continuous function of the form As(t) = —Et for 
t E [ti, £2], being constant for t G (—00, ti) and t £ (ia, +00). The effects of particle-creation 
by the T-constant field have been studied in detail in [10) . In particular, it was shown that 
in case 

T> (eE)~ 1/2 [l + m 2 /eE] (2) 

all the finite effects caused by particle-creation reach their asymptotic values, whereas the 
details involving the form of switching the field on and off can be neglected. In our calcula- 
tions, we assume this restriction from below for the time T. 

Taking the time instant t% — = T/2 — immediately before the electric field has 
been turned off, we shall now examine the mean energy value (H) of the spinor field on 
condition that its state at the initial instant t m — > — 00 should be vacuum. The potentials 
of the T-constant field do not depend on the spatial coordinates, which implies that (H) is 
proportional to the space volume V. In the one-loop approximation, (H) = wV, with the 
mean energy density w being independent of the spatial coordinates, 



w= ~(0,m| [ip(x)i ,Hip(x)] |0,m) 



(3) 



where TL is the one-particle Dirac Hamiltonian; ip(x) are the operators of the Dirac field 
in the generalized Furry representation (see, e.g., [7]) obeying the Dirac equation with the 
external background; |0,m) is the initial vacuum state in the same representation. The 
above-mentioned choice allows one to take a complete account of the pair-creation effect 
during the entire time. Also, since the electric field has not yet been switched off, this allows 
us to make a complete study of the vacuum polarization effect. Notice that the initial 
vacuum |0,m) is identical with the vacuum of those free particles that correspond to the 
initial potential A 3 = ET/2. 

The expression for w is obviously real-valued. One can see that it can be represented as 
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^ lim Q tr [(d - d' ) S in (x, x')] + ^ u ™j_ tr [(do - #o) S in (x, x')] 



(4) 



XL=x',X°=t 2 -0 



where tr [• • • ] is the trace in the space of 4 x 4 matrices, and Si n (x, x') is the so-called in — in 
Green function, 



Sin(x, x ) — i(0, in\Tip(x)tp(x )|0, in) = S c (x, x) + S p (x, x ) 



(5) 



where S c (x,x') is Feynman's causal Green function, while the function S p (x,x') is a dif- 
ference of two Green's functions, satisfying the homogeneous Dirac equation; see [2- The 
final vacuum \0,out) is the vacuum of free particles in the generalized Furry picture and 
corresponds to the constant potential A3 — —ET/2. 

The separation of Si n into the c- and p-parts is responsible for the separation of w into 



the two respective summands w 
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One can verify that the expression for w c has 



a finite limit at T — > 00, i.e., it permits a transition to the limit of a constant electric field. 
Then S c can be presented by a proper-time integral, see [11] . Using this expression, one can 
readily verify that w c is expressed in terms of the real- valued part of HEL (at B = 0). This 
contribution is due to vacuum polarization. In a superstrong electric field, it has the form 
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The contribution w p arises due to particle-creation. It is computed as follows. First of 
all, using the general theory of particle-creation (see, [7]), one can represent the function 
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S p (x,x') in the form 



%,')^E-iW [ G (+\~) G (-\~)} nm +'*Pm( X ')- (6) 
nm 

Here, {±ip n (x)} are the so-called m-solutions of the Dirac equation in a T-constant electric 
field, their asymptotics at t < t\ being stationary states of free electrons (+) and positrons 
(— ) for the Dirac Hamiltonian with the constant potential A3 — ET/2. The matrices G (±\ ± ) 
(being a matrix generalization of the Bogolyubov coefficients) are defined by decompositions 
of the so-called owi-solutions in the m-solutions: ± ip(x) — + -0(^)G( + | ± ) +_ ip(x)G (-\ ± ) . 
Here, { ± V' n (a;)} are the owi-solutions of the Dirac equation in a T-constant electric field; 
their asymptotics at t > ti describe free particles (electrons (+) and positrons (— )) with an 
energy spectrum defined by the Dirac Hamiltonian with the constant potential A3 = —ET/2. 
The matrices G (± \^) are expressed via the inner products of in- and oui-solutions, and obey 
some unitary relations following from normalization conditions for these solutions. 

In a T-constant uniform electric field, we can choose the quantum numbers of particles 
asn = (p, r), where p is the particle momentum and r = ±1 is the spin projection; then the 
matrices G (±\ ± ) are diagonal, G (±| ± )„ „/ = S r .r'dp, P 'g (±l ± ) p r ■ Here, we use the standard 
volume regularization, so that <5(p — p') — ► $ P , P '- In addition, the differential mean numbers 
of electrons (equal to the corresponding differential mean number of pairs) created from 
vacuum with a given momentum p and spin projections r are given by K P:I , = |g(_| + )| . 
It is easy to verify, using ([6]), that the function S p , which enters the expression for w p at 
x ss x' , can be taken there in the form 

S P (x,x') = -i dp ^2 ^P-r [ + ll> Ptr (x) + V> pir (x') - ~1p p , r (x) ~^p, r { x ')] > x ~ x ' ■ 
J r=±l 

For t > t2, the oi/i-solutions ^i^ p r (x) describe free particles with the quantum numbers p, r 

and energies e p , r = J m 2 +p\ + (eET/2 — p 3 ) 2 , p^ = (p^p^O). Taking all the above 
into account, we obtain 

wP = 4^3 / d P E ^S Pir ■ (7) 

« LI 



r=±l 



This quantity is the mean energy density of pairs created from vacuum. It can be estimated 
in the case of strong electric fields, E > E c and a sufficiently large T as follows. As has been 
demonstrated in [10], in case the time T is sufficiently large, 

T > (m 2 + pi + eE) {eE)^' 2 , 

and the longitudinal momenta are restricted by the condition |_p3 1 < (V eET/2 — K^j V eE, 

2 ' 



where K p is a sufficiently large arbitrary constant, v eET ^> K p ^s> 1 + (m + pjj /eE, the 
differential mean numbers H PjT . have the form 



Np,r = exp 



m 2 + p\ 

eE 



For any fixed pj_, the function H Pir is fast-decreasing for \p%\ > eET/2 — K p j \J eE. For 

this reason, we can disregard the contribution to the integral ([7]) due to the integration 
over such momenta p% in comparison with the main contribution, which is defined by the 
dimcnsionless parameter eET 2 . The latter parameter, in fact, determines a large integration 
domain over p%. In its turn, the exponential decrease of Kp^ with the grows of p^ allows one 
to ignore the contributions to the integral (|7|) due to a large p\/eE > \fe~ET . Consequently, 
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in order to evaluate the term which leads in y/eET in integral |(7J) we can replace H Pjr 
by Np S ^ under condition ([2]), while restricting the limits of integration over momenta by 

the region \p 3 \ < VeE (VlET/2 - K^j , where K is a sufficiently large arbitrary constant, 

sfeET K ^> 1 + rn 2 /eE . Having calculated the integral (JT]) over p±, we obtain the 
T-leading term in the form 

vT^eETX, H^—^exp^r— j . ( 8 ) 

III. We now suppose that the energy density of particles w — w p , arising precisely due 
to the action of a T-constant electric field, should be essentially smaller than the density of 
the electric field itself, being equal to the classical Maxwell density of energy, £^°> = E 2 /8n. 
Thus, the condition of a smallness of back-reaction is w p <C E 2 /8tt, which, owing to ([5]), 
takes the form of a restriction from above on the dimensionless parameter eET 2 : 



e^«-ex P ^-j. (9) 

On the other hand, all the asymptotic formulas have been obtained under condition @, 
which restricts the mentioned parameter from below, [l + m 2 /eE] < eET 2 . Since 7r 2 /2e 2 ^> 
1, there exists a region of values of E and T that satisfies both the inequalities. We note 
that time scale from above in ^ is more restrictive than the scale derived from the rate of 
pair production, see [12] . 

In case the initial state is in thermal equilibrium at temperature 9, the mean energy 
density has an additional term Wg, which represents, in fact, the work of a T-constant field 
on particles from the initial state, as well as the term Wg 

Wg = -7-3 [ d P ^2 ^p.r n P,r (in) £ p , r , n p , r [in) = [exp (e Ptr /6) + , 

n J r=±l 



where e Pir = y m 2 + + {qET/2+p 3 y is the energy of a free m-particle. The latter 
term determines a temperature-dependent correction to the energy of particles created from 
vacuum; see [9]. The energies of particles that contribute to Wg in the limit of a large T 
are mostly determined by a large longitudinal kinetic momentum, with the energy being of 
order eET, as well as the energies of particles created at 8 — in the expression ([5]) for w p . 
Given that, however, the density of initial particles is constant, being determined only by 
the initial condition, whereas the density of created particles increases in proportion with T. 
Therefore, at large T and E, Wg can be neglected in comparison with w p , and w ~ w p + Wg . 

In case the initial state is in thermal equilibrium at low temperatures 9 <C eET, the 
contribution Wg turns out to be small in comparison with w p . At high temperatures 9 3> 
eET, the energy density has the form w — (eET/69) w p . Thus, the restriction ([9]) is valid 
both for the vacuum initial state and for a low-temperature initial thermal state. At high 
temperatures we have a weaker restriction: 



(eE) 2 T 3 ^ 3tt 2 

I? ~" r V " eE 



<— exp(7r— ). (10) 



Analogously, one can find restrictions for QED with charged bosons in a T-constant 
electric field. At low temperatures, we have 

7T 2 f m 2 \ 
e£T 2 « - exp ( , 
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where J is the number of the spin degrees of freedom (J = 1 for scalar particles and J = 3 for 
vector particles). In the case of high temperatures the restriction has a completely different 
character than (fit)!) , namely, 

.Tln(^T)«^ex P ( 7 r^). 

One can easily extend these results to D + 1 dimensions, using the corresponding AT in ((SJ), 
taken from Eq. (37) in [ID] . 

IV. A similar analysis can be performed in the case of QCD with an electric-like non- 
Abelian external background. Such a background is a part of the known chromoelectric 
flux-tube model [13] ■ At present, the chromoelectric field is associated [2] with an effective 
theory, color glass condensate. Here, we shall derive restrictions on the external background 
which allows one to treat particles created from vacuum still as weakly coupled, owing to the 
property of asymptotic freedom in QCD. To this end, we use the results obtained in |15[ll6j 
for QCD with a constant ST/ (3) chromoelectric field E a (a = 1, ... ,8). If the initial state 
is vacuum, the density of created gluons is noticeably higher than the density of created 
quarks at any intensity of a T-constant chromoelectric field; see [15]. The same is valid at 
any finite temperature; therefore, for our purposes it is sufficient to take into account only 
the gluon contribution. It has been demonstrated in [15] that the pj_-distribution density 
n giuon Q £ g mong produced from vacuum with all the possible values ps and the quantum 
numbers that characterize the inner degrees of freedom can be presented as follows: 

<T n = ^jy^)^, N«) = exp - (11) 

where E/j-s are positive eigenvalues of the matrix if abc E c for the adjoint representation of 
SU(3), and q is the coupling constant. The (j)-terms in dTTI) can be interpreted as those 
obtained for Abelian-like electric fields E^, respectively. Then, the total energy density 
of gluons created from vacuum by the field E^ is determined by integrals of the kind ([7]). 
Taking into account that maxima of the fields are restricted by the condition Er^ < \[C\ 
(C\ = E a E a is a Casimir invariant for SU(3)) and the relation Y^=i ^(j) ~ 3Ci/2, one can 
find that at low temperatures 9 <C q\[C\T the consistency restriction for the dimensionless 
parameter q\fC\T 2 has the form 

qyfGxT 2 < ir 2 /3q 2 . 

As in the case of QED, this restriction must be accompanied by a restriction from below, 
1 <C q^/CiT 2 , which is related to the fact that all the asymptotic expressions have been 
obtained for sufficiently large values of T. Therefore, the T-constant SU (3) chromoelectric 
field approximation is consistent during the period when the produced partons can be treated 
as weakly coupled, due to the property of asymptotic freedom in QCD. At high temperatures, 
9 3> q\fC\T, the consistency restriction is far more rigid: 

9T\n (qV~C~iT 2 ) « 7r 2 /3g 2 . 

The above established consistency restrictions determine, in fact, the time scales from 
above of depletion of an electric field due to the backreaction. 
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